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ABSTRACT 

New HST WFPC2 observations of the lensed double QSO 0957+561 will allow improved constraints 
on the lens mass distribution and hence will improve the derived value of Ho. We first present improved 
optical positions and photometry for the known components of this lens. The optical separation between 
the A and B quasar images agrees with VLBI data at the 10 mas level, and the optical center of the 
primary lensing galaxy Gl coincides with the VLBI source G' to within 10 mas. The best previous model 
for this lens (Grogin and Narayan 1996) is excluded by these data and must be reevaluated. 

Several new resolved features are found within 10" of Gl, including an apparent fold arc with two 
bright knots. Several other small galaxies are detected, including two which may be multiple images of 
each other. We present positions and crude photometry of these objects. 

Subject headings: distance scale — gravitational lensing 



1. INTRODUCTION 

Now that an accurate time delay between the two im- 
ages of the gravitationally lensed quasar Q0957+561 has 
been reliably determined (Kundic et al. 1996), it should be 
possible to use this system to measure Ho to an accuracy 
of a few percent (Refsdal 1964). This measurement would 
be independent of the usual distance ladder and insensi- 
tive to local motions. The largest remaining uncertainty 
in this method is our ignorance of the mass distribution 
in the lens. This ignorance takes two forms: the first is 
a degeneracy in the lens models under the addition of a 
smooth sheet of matter across the central 10-20" of the 
lens. Changing this smooth component in a mass model 
leaves all image positions and magnifications unchanged, 
yet alters the derived Ho. Our first step in reducing the 
uncertainty in the value of Ho was to obtain a weak lens- 
ing estimate of the smooth mass distribution contributed 
by the small galaxy cluster at the primary lens redshift 
(Fischer et al. 1996). Recent measurements of the veloc- 
ity dispersion of Gl may also help reduce this degeneracy 
(Falco et al. 1997). 

Our second form of ignorance is of the structure of the 
mass distribution in the strong lensing region within 10" of 
Gl. Weak lensing statistical mass reconstructions have in- 
sufficient resolution to determine structure on these small 
scales. Existing models of the lens presume some para- 
metric form for the mass distribution of the primary lens 
galaxy Gl, and adjust the parameters until the quasar 
image positions are reproduced {e.g. Falco, Gorenstein, 
& Shapiro 1991; Bernstein, Tyson, & Kochanek 1994 
[BTK94] ; Grogin & Narayan 1996). Tighter and more 



numerous constraints on the optics of the lens will lead to 
more detailed and precise models of the mass distribution, 
and hence more accurate values of Ho . The best existing 
constraints are the highly accurate quasar core positions 
from VLBI observations (Gorenstein et al. 1988) and the 
further VLBI observations of the structure of the 50 mas 
jets which extend from each quasar core (Garrett et al. 
1994). We have obtained deep WFPC2 observations of 
the 0957+561 system to see if the visible-light positions of 
the quasars match the VLBI positions, and to see which 
(if any) of the reported weak radio sources corresponds to 
the center of Gl. Furthermore, nature has been kind, and 
the WFPC2 images reveal other objects in the background 
of Gl which are strongly lensed, allowing us to place new 
constraints on the mass distribution. 

2. OBSERVATIONS 

Observations of the field around the double QSO 
0957+561 were taken with the WFPC2 on 19 Nov. 1995 
and 26 Nov. 1995 with the QSO centered in WF3. Two fil- 
ters were employed, with 2 x 80s and 14 x 2300s in F555W 
and 2 x 160s, 1 x 500s and 2 x 900s in F814W. The images 
were pipeline processed in the standard manner. The ob- 
servations were dithered by an integer number of WF pix- 
els with a maximum offset of ±2" in either dimension. The 
fourteen 2300s exposures were shifted, combined and sky 
normalized using the Tukey biweight algorithm (Hoaglin, 
Mosteller, & Tukey 1983), which eliminates cosmic rays 
on the final combined image. Figure la is the shallow 
80 sec F555W image. Figures lb and lc are the deeper 
combined F555W image and are discussed below. 
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Fig. 1. — Panel (a) is the 0957+561 lens system as seen in 80 second F555W WFPC2 images, with the A and B quasar 
images and primary lensing galaxy apparent. Panel (b) shows a combine 39 ksec F555W image; a shifted version of the 
image has been subtracted so that the A quasar acts as PSF template for B quasar. The newly discovered arc and "blobs" 
are marked. White objects are negative "ghosts" from the subtraction process. Orientation and scale are same as panel 
(a). Panel (c) shows the same image after subtraction of the fitted elliptical isophotcs for Gl. The arc and its bright spots 
arc more apparent; other features in the vicinity of Gl are residuals from the B and Gl subtractions. The inset shows 
zoomed images of Blobs 2 and 3, which may be multiple images of the same background galaxy. 
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3. ASTROMETRY 

With the resolution of the HST we should be able to de- 
termine whether the optical quasar positions coincide with 
the cores of the VLBI radio sources. More importantly for 
lensing models, we wish to determine the position of the 
center of galaxy Gl, and sec whether it coincides with ei- 
ther the VLBI object G' (Gorenstein et al. 1988) or the 
VLA object G (Roberts et al. 1985). 

3.1. Quasar Separation 

The centroids of the A and B quasars were determined 
on each of the four short exposures, and the STSDAS 
task METRIC was used to remove the geometric distor- 
tions on the WF3 chip and place the objects on the J2000 
system. The vector from the B to the A quasar image 
has measured length rnsT — 6'.'169(3), at position angle 
PA H st = -H?4(l) (from N through E). The uncertain- 
ties are dominated by calibration factors: the pixel scale 
of WF3 is uncertain by « 3 parts in 10 4 (Cox 1996), and 
the PA is uncertain because the HST roll angle is calcu- 
lated from the positions of the two guide stars, which are 
known only to « 0.5" accuracy. 

The most accurate measures of the (A-B) separation are 
from the VLBI data of Gorenstein et al. (1988), who deter- 
mine r VLB i = 6'.'17499(2) and PA VLB i = -1T?7029(2) in 
epoch 1950. After rotating the PA into J2000, the VLBI 
data yield r VLBI = 6'.'175 and PA VLBI = -11?453, with 
miniscule error compared to the HST data. The HST and 
VLBI distances from B to A agree to 6 mas (2a), one part 
in 10 3 . The PA of the (A-B) vector agrees to within our 
0?1 knowledge of the HST roll angle, which corresponds 
to a 10 mas uncertainty in the (A-B) vector. Because 
the agreement is good and the uncertainties in the optical 
positions are dominated by calibration factors, we will as- 
sume henceforth that the optical positions coincide with 
the VLBI cores, and will use the VLBI (A-B) vector to 
calibrate the orientation and scale of the HST image. The 
uncertainties in the HST guide star catalog preclude a use- 
ful comparison of the absolute positions of the quasars in 
the optical vs. radio images. 

3.2. Gl Position 

The lensing galaxy Gl is visible on each short HST ex- 
posure, but the low S/N and preponderance of cosmic rays 
make it more profitable to determine its center from the 
combined deep HST image. To successfully measure Gl, 
it is necessary to subtract the B quasar flux from the im- 
age. We find that the A quasar image makes a far better 
PSF template than can be produced by any other method. 
From the combined deep image we subtract a copy which 
is shifted so that the A quasar overlapps B. yielding the 
image shown in Figure lb. 

A robust estimate of the Gl center was made by cen- 
troiding and ellipse-fitting over different apertures on the 
combined image. We cannot measure the positions of A 
and B on the combined deep image because they are sat- 
urated, so we must transfer the short-exposure coordinate 
system (as determined by forcing the (A-B) vector to agree 
with VLBI data) onto the combined image using the 8-10 
objects which are visible on the short exposures yet un- 
saturated in the deep images. The errors in this transfer 



are comparable to the errors in the Gl center. The optical 
position of the center of Gl is (0"1776(35) E, l'/0186(35) 
N) (J1950) of quasar B. This is fully consistent with the 
ground-based position for Gl (0.19 ±0.03" E, 1.00 ±0.03" 
N) of B (Stockton 1980). We list J1950 positions in this 
section for easy comparison with previous works; the J2000 
positions of the various central objects are listed in Table 1. 

The center of Gl is within 10 mas (2.7a) of the VLBI 
source G' at (G'-B) = (0"181(l) E, l'.'029(l) N) reported 
by Gorenstein et al. (1983). Independent EVN observa- 
tions at 18 cm by Garrett (1990) place G' at (O 7 / 179(1) E, 
1"026(1) N) relative to B (all are J1950). Thus it seems 
that the G' radio source is most likely at the center of Gl; 
the 10 mas difference between G' and Gl centers, if real, 
is only ss 40/i _1 pc at the distance of Gl. If G' were the 
third image of the quasar core, we would expect it to be 
demagnified relative to the A and B images by a factor 
~ rQi/rE, where re is 10 mas and r B is the Einstein ra- 
dius of the lens, roughly 3". Thus G', if a third image, 
should have a flux ^0.3% that of the B image, whereas 
Garrett (1990) measures a G'/B flux ratio of 2.8%. Similar 
G'/B flux ratios were found by Gorenstein et al. (1983), 
suggesting that G' is about an order of magnitude too 
bright to be the third image of the quasar. 

The radio source G found in 6 cm VLA maps of the sys- 
tem (Roberts et al. 1985) is significantly displaced from 
both our optical and the VLBI sources: (G-B)=(0"151(l) 
E, 1"051(1) N). Likewise 18 cm MERLIN observations 
by Garrett (1990) yield (G-B)=(0'.'148(8) E, 1"093(8) N). 
Both of these observations find G many sigma to the North 
and West of G' and Gl. The discrepancy is likely due 
to structure in the radio sources which is not resolved by 
the VLA or MERLIN observations; the VLBI observations 
clearly show such structure in the A and B images, and 
there may be jets in G as well which are undetected by 
VLBI. These issues are discussed in detail in the radio 
observation papers. 

Finally we note that the lens model for 0957±561 by 
Grogin & Narayan (1996) places the lens (Gl) center at 
(0"215 E, 1"057 N) of B, 15cr from the new optical position. 
These authors note that their model fits the quasar jet 
positions extremely poorly if the galaxy is forced to reside 
near G', as indicated by the HST images, so a reevaluation 
of the lens models is now required. 

4. SURFACE PHOTOMETRY OF Gl 

The surface photometry task ISOPHOTE in the STS- 
DAS package is used to fit elliptical isophotes to Gl. 
Isophote fitting is done on the image shown in Figure lb, in 
which the A quasar has been shifted and subtracted from 
the B quasar. Note that this leaves a negative "ghost" 
quasar 6" south of Gl. This ghost and the A quasar im- 
age limit the radius to which surface photometry can be 
done. During the ellipse fit we mask the areas covered by 
the B quasar core and the newly discovered small objects 
described in §|| Figure lc shows the WFPC2 image after 
subtraction of the fitted model for Gl. 

The surface brightness, ellipticity, and PA of the Gl 
isophotes are plotted vs radius in Figure 2. We take 
the measured F555W surface brightnesses to be equiva- 
lent to V band, since Fukugita et al. (1995) calculate 
F555W-V for a z = 0.36 elliptical to be < 0.03 mag. We 
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have overplotted the ground-based i?-band surface pho- 
tometry of Gl from BTK94, assuming the predicted color 
V — R = 1.55 that a present-day elliptical galaxy would 
have if observed at z = 0.36 (Coleman, Wu, and Weedman 
1980). This leaves the BTK94 data « 0.2 mag fainter than 
the WFPC2 data in the region of overlap. This small dis- 
crepancy may be ascribed to passive evolution in Gl and 
perhaps calibration errors. 

The WFPC2 and KPNO datasets match smoothly in all 
measured parameters. The HST surface photometry in- 
side 2" radius continues the trends seen farther out in the 



BTK94 profiles — a slight isophotal twist and a tendency 
toward rounder isophotes in the center of Gl [wiggles in 
the ellipticity and PA profiles are probably due to residual 
quasar flux]. The surface brightness profile is no longer 
well fit by a single power law or by a deVaucoleurs pro- 
file; the latter is too flat in the center of the galaxy. Like 
most ellipticals, Gl continues a power-law rise in surface 
brightness to the center at the resolution limit of the HST 
(Gebhardt et al. 1996), which means that the third quasar 
image is likely to be highly demagnified. 




Semi-Major Axis (arcsec) 



Fig. 2. — Fitted elliptical isophotes for Gl are plotted vs the semi-major axis of the isophotes. Filled symbols are from 
fits to the WFPC2 image F555W image, which should be very close to V^-band. Open symbols are from the ground- 
based i?-band image described in BTK94. The upper panel shows the surface brightness profile; the R data are plotted 
assuming V — R = 1.55 as expected from nearby elliptical galaxies redshiftcd to z = 0.36. Middle and lower panels show 
the ellipticity and PA of the fitted isophotes; most of the scatter is due to residuals from the subtraction of quasar B, but 
in general we see good agreement between the 2 data sets, a trend toward rounder isophotes in the center, and a slight 
isophotal twist. The "A" and "B" at the top axis mark the radii of the A and B quasar images. 



5 



5. NEWLY DISCOVERED FEATURES 

Several previously undiscovered objects are visible 
within 10" of Gl on the quasar-subtracted WFPC2 image 
(Figure lb) and the galaxy-subtracted image (Figure lc). 
At least some of these are likely to be lensed background 
objects, which will place new constraints on lens models 
for this system. Table 2 lists positions, total magnitudes, 
and peak surface brightnesses (at the WFPC2 resolution) 
of the new features. 

5.1. Arc 

A previously undetected "arc" is visible just outside the 
B QSO image. The arc is visible even before the subtrac- 
tion of the B QSO quasar flux, thus is not an artifact of 
the PSF subtraction. There appears to be substructure in 
the arc, specifically two bright spots ("knots") separated 
by 0'.'5. Rough coordinates of points along the arc and of 
Knots 1 & 2 are given in Table 2. The arc is of quite low 
surface brightness and the knots are faint (V w 26) and 
perhaps unresolved. Because of the low S/N of these fea- 
tures and the difficult QSO and galaxy subtraction, the 
uncertainties on the knot positions and fluxes are only 
crudely estimated. 

Unlike the putative arc shown in BTK94, this candi- 
date arc looks very much as one would expect for a fold 
arc given the geometry of the QSO images. Arcs are found 
astride critical lines, and for sources at redshifts well be- 
hind the z = 0.36 lens there must be a critical line at a 
radius from Gl intermediate to the radii of the A and B 
quasar images. Indeed the new arc candidate is found at 
such an intermediate radius. The parts of a fold arc inte- 
rior and exterior to the critical line are mirror images of 
the same section of the source. It seems likely that Knot 
1 and Knot 2 are multiple images of the nucleus of or a 
bright spot in a faint source galaxy. The two knots have 
similar fluxes and similar peak surface brightnesses, con- 
sistent with a fold arc scenario. Higher-S/N images will 
show more detail in the arc and knots, permitting a more 
definitive conclusion as to the location of the critical line. 

5.2. Small Sources 

Seven other sources are visible in the vicinity of the A 
and B quasar images. These "blobs" are labelled in Fig- 
ure lb and their positions, fluxes, and peak surface bright- 
nesses are are listed in Table 2. Uncertainties are difficult 
to quantify because these are mostly extended objects and 
the sky subtraction is the leading source of error. We 
list all positional uncertainties as 0'.'05, or one-half of a 
WFPC2 pixel. 

Blob 1 is unresolved. It is likely to be at the Gl red- 
shift or a foreground star, because if it were behind Gl 
it would almost certainly be multiply imaged. No other 
stellar objects are visible in the strong lensing region. 

Blob 2 is only l'/4 from Gl and hence will be multiply 
imaged if at redshift higher than Gl. Blob 2 is clearly 



resolved and hence non-stellar. Since the median redshift 
for V = 25 galaxies is well above z = 0.35, it is very likely 
that a second image of Blob 2 is present in the field. A lens 
model which we fit to the quasar positions suggests that 
Blob 3 has the correct position and flux to be this second 
image. The peak pixel of Blob 3 has ~ 2 times higher 
flux than the peak pixel in Blob 2, but this difference in 
apparent surface brightness could be due to an unresolved 
bright spot in the source galaxy. As can be seen in the 
inset to Figure lc, both galaxies are extended in the EW 
direction and are brighter at their W ends; the lens model 
suggests that the counterimage of Blob 2 should be flipped 
about a nearly horizontal axis, so their crude morphologies 
are consistent with multiple imaging. Higher-S/N imag- 
ing will again allow a more definitive test of the multiple 
imaging hypothesis. 

The other four new objects are sufficiently distant from 
Gl that they are probably not multiply imaged regardless 
of their redshifts. 

6. CONCLUSIONS 

The WFPC2 images provide a wealth of informa- 
tion useful in modelling the strong-lensing effects in the 
0957+561 system. The optical A and B quasar and Gl 
galaxy positions agree with the enormously precise VLBI 
positions for A, B, and G' to within 10 mas. This forces a 
revision of the best published model of the 0957+561 lens 
(Grogin and Narayan 1996). Fortunately the WFPC2 im- 
ages provide new constraints to use in modelling the sys- 
tem: the arc constrains the location of the critical line, es- 
pecially if we can interpret the two bright knots as multiple 
images; and the hypothesis that Blobs 2 and 3 are multiple 
images of a common source adds another set of deflection 
and magnification constraints to candidate models. These 
new data add at least 4 new constraint equations to mod- 
els of this lens. In a later publication we will explore lens 
models which can satisfy these new constraints as well as 
the position and flux constraints from previous radio imag- 
ing. The new STIS CCD imager on HST will be trained 
on this system in Cycle 7 to provide higher S /N and better 
PSF sampling on the newly discovered arc and background 
galaxies. These images should show sufficient detail in the 
arc and faint background galaxies to definitively test the 
lensing hypothesis, to determine the full relative magnifi- 
cation matrix for each lensed pair, and perhaps to detect 
the fainter counter-images of the arc, further increasing 
the constraints on the model. This should allow exquisite 
accuracy in the determination of the mass distribution in 
this system, and improve the accuracy of the determina- 
tion of Hq. 

Support for this work was provided by NASA through 
grants # HF-01069.01-94A (PF) and #GO-05979.0X-94A 
(GB, GR, JAT) from the Space Telescope Science Insti- 
tute, which is operated by the Association of Universi- 
ties for Research in Astronomy Inc., under NASA contract 
NAS5-26555. 
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Table 1 

Positions for Central Object 



Object 


Instrument 


RA 1 


Dec 1 


Uncertainty 


Reference 


Gl 


HST 


(X'1820 


l'/0178 


0"0035 


This work 


Gl 


CFHT 


Of! 19 


l'/OO 


Of! 03 


Stockton 1980 


G' 


VLBI 


Of! 185 


l'/028 


07001 


Gorenstein et al. 1983 


G' 


EVN 


07183 


l'/025 


O'.'OOl 


Garrett 1990 


G 


VLA 


Of! 155 


l'/050 


O'.'OOl 


Roberts et al. 1985 



displacement from B center (J2000) listed. 



Table 2 

Faint Object Positions and Photometry 



Object 


RA 1 


Dec 1 


Magnitude 2 


Peak SB 2 


Remarks 


Quasar A 


-l'/408 


+5'.'034 








Quasar B 


-07182 


-T/018 








Arc 


-l'/51 


-l'/81 




;>24.7 


Point along arc 


Arc 


-or 93 


-2'.' 21 




;>24.7 


Point along arc 


Knot 2 


-0'/48 


-2"43 


25.8 ±0.5 


23.4 ±0.1 


Position of peak flux 


Knot 1 


-0'/06 


-2'! 55 


26.1 ± 0.5 


23.6 ±0.1 


Position of peak flux 


Arc 


+0'/47 


-2'/66 




;>24.7 


Point along arc 


Arc 


+T/09 


-2'.' 74 




;>24.7 


Point along arc 


Blob 1 


-T/95 


-0'/46 


25.6 ±0.1 


22.26 ±0.04 


Unresolved 


Blob 2 


+l'/54 


-0'/05 


24.8 ±0.1 


22.95 ±0.07 


Possible counterimage of Blob 3 


Blob 3 


-2'.'86 


+3'/47 


23.5 ±0.1 


22.32 ±0.04 


Possible counterimage of Blob 2 


Blob 4 


-2'/70 


+6'/54 


25.6 ±0.1 


22.64 ±0.05 




Blob 5 


-5'/67 


+4'/58 


26.4 ±0.1 


23.6 ±0.1 




Blob 6 


-7'/22 


+2'/77 


25.6 ±0.1 


23.2 ±0.1 




Blob 7 


+3'/65 


+5'/66 


24.8 ±0.1 


22.48 ±0.04 





displacement from Gl center (J2000) listed. All positions except Quasars uncertain by ±0'.'05. 
2 F555W magnitudes are given; V mag should differ by < 0.05 mag. 



